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1 . INTRODUCTION 


The  singularity  expansion  method  (SEM)  has  been  used  to  express  the 
transient  response  of  antennas  and  scatterers  in  a compact  and  useful 
form.*  By  usinq  Prony's  algorithm,  Van  Blaricum  and  Mittra1  were  able 
to  extract  the  SEN  poles  and  residues  from  analytical  waveforms.  The 
singularity  expansions  of  experimental  data  have  also  been  obtained,* 
and  Van  Blaricum  and  Schaubert'  have  demonstrated  that  the  SEM  can  be 
used  effectively  to  compute  an  antenna's  impulse  response  from  measured 
data . 

The  use  of  the  SEM  to  calculate  equivalent  impedance  circuits  has 
also  been  discussed  by  many  authors.^  However,  there  has  been  no 
previous  application  of  the  SEM  to  calculate  equivalent  lumped-element 
impedance  circuits  for  antennas.  Such  circuits,  when  calculated  by 
usinq  the  SEM,  will  be  valid  over  a wide  frequency  range  and  will, 
therefore,  be  extremely  useful  for  transient  and  out-of-band  response 
analyses.  Additional  advantages  of  the  SEM  approach  are  that  the 
goodness  of  fit  to  the  data  can  be  judged  prior  to  calculating  the 
circuit  element  values,  and  the  impedance  function  can  be  readily 
synthesized  by  existing  techniques. 

Several  problems  relevant  to  the  computation  of  equivalent  circuits 
for  antennas  have  been  investigated,  and  the  results  are  presented  in 
sections  2,  3,  and  4.  Throughout  the  investigation,  primary  importance 
has  been  placed  upon  utilization  of  experimentally  or  numerically 
derived  data.  The  SEM  expansion  of  time-domain  ref lactometer  (TDR)  data 
for  a dipole  antenna  with  a quarter-wavelength  ( V4)  balun  has  been 
obtained  and  used  to  calculate  a lumped-element  equivalent  circuit. 


2.  RELATIONSHIP  BETWEEN  Z(s)  AND  MEASURED  QUANTITIES 

The  first  step  in  calculating  an  equivalent  circuit  is  deriving  the 
relationship  between  the  antenna  impedance  function  Z(s) , where  s is  the 
Laplace  transform  frequency  variable,  and  a measurable  quantity.  One  of 
the  most  common  and  most  easily  used  instruments  for  measuring  impedance 
is  the  TDR,  which  excites  the  antenna  with  a step  voltage  and  displays 

L.  Van  Blaricum  and  R.  Mittra,  A Technique  for  Extracting  the 
Poles  and  Residues  of  System  Directly  from  Its  Transient  Response, 
IEEE  Trans.  Antennas  Profvtg.,  AP-23  (November- 1975) , 777-762. 

•M.  L.  Van  Blaricum  and  P.  H.  Schaubert , An  Exf>erimental  Transient 
Transfer  Function  via  Prony's  Method,  Proc.  1976  VSNC/VRSI  Meeting, 
Amherst,  MA  (October  1976),  97. 

*See  Selected  Bibl iography — Transient  Response  of  Antennas  and 
Scatterers . 

t See  Selected  Bibl iography — Experimental  Pata. 

tSee  Selected  Ribl iography — Equivalent  Circuit  Synthesis. 
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the  resultinq  terminal  voltage  as  a function  of  time  (fig.  1) . This 
terminal  voltage,  v(t),  is  related  to  the  antenna's  impedance  through 
the  reflection  coefficient  and  a convolution  integral.  Specifically,  if 
the  Laplace  transform  iv»ir  is  defined  by 


■ st 


’(a)  - | v(t)e  dt, 

•'o 

,(t)  * ds  • 


(la) 


(lb) 


where  a is  an  appropriate  constant,  and  j - *-l,  then 


V (s)  - -v  (s)  (1  ♦ R(s)  j . 

2 0 


(2) 
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Figure  1.  Time-domain  ref lectometer  measurement 
of  antenna  impedance. 

In  equation  (2),  V(s)  and  V^(s)  are  the  Laplace  transforms  of  the 
terminal  voltage  and  excitation  voltage,  respectively,  and  R(s)  is  the 
voltage  reflection  coefficient  of  the  antenna, 


Z(s)  - Z, 


R(s)  - 


Z(s)  + Z 


(3) 


0 


where  Z is  the  source  impedance  of  the  TPR. 


Applying  the  inverse  Laplace  transform  (lb)  to  equation  (2)  yields 


j 


v(t)  - - / v (t  - t ) [ 6 < t ) + R ( t ) ) dt  . (4) 

A A)  W 

Since  R(t)  is  the  inverse  Laplace  transform  of  equation  (3),  the 
relationship  between  v(t)  and  Z(s)  is  not  trivial.  In  the  Laplace 
frequency  domain,  however,  the  relationship  is  quite  simple. 
Substituting  equation  (3)  into  equation  (2)  yields 


V(s) 


V0<S>2 


Z(s) 
(s)  t 


Z0 


That  is, 


Z(s) 


7 V(3> 

0 V0(s)  - V (s)  • 


(5) 


(6) 


It  is  clear  that  an  analytical  expression  for  Z(s)  can  be  obtained  if 
analytical  expressions  for  V,(s)  and  V(s),  which  are  the  Laplace 
transforms  of  measurable  quantities,  are  available.  The  SEM,  together 
with  Prony's  method,  provides  an  efficient  means  of  obtaining  analytical 
expressions  for  V,(s)  and  V(s)  when  v0(t)  and  v(t)  arc  known  in  sampled 
data  form. 


3.  SINGULARITY  EXPANSION  AND  PRONY'S  METHOD 

The  singularity  expansion  of  a transient  waveform  v(t)  is  defined  as 
the  exponential  approximat ion 


v<t>  V n • (7) 

n=l 

This  expansion  can  be  readily  transformed  to  yield 


V (s) 


<a> 


The  s are  the  poles  of  the  expansion  and  the  A are  the  residues. 
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TTie  problem  of  importance  in  this  section  is,  given  v(t)  at  discrete 
points  t « t^,  i = 1,  2,  . . . , find  the  A and  s that  provide  the  best 
approximation  in  the  form  of  equation  (7).  Because  the  s are  not 
known,  this  is  a nonlinear  approximation  problem.  However,  nProny  has 
developed  an  algorithm  that  converts  this  problem  into  two  linear 
problems  plus  one  polynomial  root  finding  problem,  which  can  be  readily 
solved  on  a digital  computer. 

An  outline  of  Prony's  method  is  presented  below.  Readers  unfamiliar 
with  the  method  may  consult  one  of  the  references.1.* 

If  the  waveform  v(t)  is  sampled  at  uniform  intervals  At,  then 
equation  (7)  becomes 


s At 

Defining  Z^  = e , equation  (9)  becomes 

N 

Vk  * £ Vn  • * - 2 <10, 

n=  1 

It  can  be  shown  that  Z satisfying  equation  (10)  are  the  roots  of  the 
polynomial  equation  n 

a0  + “l2  + a2Z2  + • • . + Zn  = 0 , (11) 


1M.  L.  Van  Blaricum  and  R.  Hittra,  A Technique  for  Extracting  the 
Poles  and  Residues  of  a System  Directly  from  Its  Transient  Response , 
IEEE  Trans.  Antennas  Propag.,  AP-23  (November  1975),  777-782. 

*See  Selected  Bibliography — Transient  Response  of  Antennas  and 
Scatterers . 
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where  the  ^ are  the  solutions  of  the  linear  set  of  equations 


(12) 


Prony's  method,  therefore,  consists  of  three  steps: 

a.  Solve  the  linear  set  of  equations  (12)  for  the  coefficients  . 

b.  Obtain  the  N roots  of  the  polynomial  equation  (11)  and, 

therefore,  the  s « (In  7.  )/At. 

n n 

c.  Solve  the  linear  set  of  equations  (10)  for  the  residues  A (7. 
are  known  from  step  b) . 

Since  v(t)  is  real,  the  poles  s (and  the  roots  7 ) must  either  he  real 
or  occur  in  complex  conluqate  pairs. 

4.  UNIFORM  NORM  APPROX  I MAT  I ON 

Once  the  poles  s are  known,  step  c of  Prony's  method  is  a linear 

approximat ion  problem.  That  is,  one  must  find  the  A that  provide  the 

best  approximat ion  to  v(t)  for  a qiven  set  of  s . Most  workers  in 

n 

electromagnet ic  scattering  have  used  the  L2  (least  squared  error)  norm 
when  calculating  the  A . This  norm  has  the  advantage  that  it  leads  to 
a linear  set  of  algebraic  equations  for  the  residues  A . However,  the 
Lt  norm  has  the  disadvantage  that  it  is  not  sensitive  to  the  large 
deviations  that  often  occur  near  t • 0.  Figure  2 shows  a singularity 
expansion  obtained  by  the  I,,  norm  approximat  ion  method. 
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Figure  2.  Typical  l^-norm  approximation  of 

time-domain  ref lectometer  (TDR)  data. 

In  order  to  improve  the  accuracy  of  calculated  singularity 
expansions,  the  third  step  of  Prony's  method  has  been  replaced  by  the 
uniform  norm  approximation  problem:  Given  the  poles  s , determine  the 
residues  A^  that  minimize  ||e  (t)  ||  where  n 

N s t 

||e  (t)  ||  = m^X  |v(t)  - £ A e n | . (13) 

n=l 

The  norm  ||  • ||  defined  by  equation  (13)  is  the  uniform  norm  and  is 
sensitive  to  any  large  deviations  of  the  approximation.  However, 
computation  of  the  best  uniform  norm  approximation  is  not  a linear 
problem.  Therefore,  additional  computation  time  is  required.  A 
comparison  of  the  L2  and  (uniform  norm)  approximation  problems  for 
ttep  c of  Prony's  method  is  contained  in  table  I.  The  two  waveforms 
that  were  used  for  this  comparison  are  shown  in  figure  3,  and  the 
parameters  are  defined  by  the  discrete  problem 

. N 

min  min  max  __  s kAt 

An  ||  e(kAt)  11^=  A^  0 <k  <.  K | v(kAt)  - £ A^  n | . (14) 

n=l 
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TABLE  I.  COMPARISON  OF  L?  AND  L„  NORM  APPROXIMATIONS 


Number  of 
points,  K 


Number  of 
poles,  N 


Computation  time  (s) 


Squared  error 


Type  of  data 


Impulse  response 

Impulse  response 

Impulse  response 

TOR 

TDR 

TDR 

TDR 


8 12 
TME  (ns) 


1 2 3 4 5 

TIME  (ns) 

(b) 

Test  waveforms  for 
Loo- norm  approximate 
response  of  dipole 
and  (b)  time-domain 
dipole. 


Figure  3 


- — rrrr 


Although  the  computation  time  is  significantly  greater  for  the  L^-norm 
approximation  than  for  the  L;-norm  approximation,  the  errors  are 
sufficiently  less  to  make  the  L -norm  attractive.  In  fact,  ths 

kU 

uniform  norm  often  yields  a very  good  approximation  when  the  Lg  norm 
totally  fails  due  to  ill-conditioning  of  the  matrix  equation. 

Since  the  uniform  norm  is  equivalent  to  the  1^,  norm  whan  p ®, 
the  above  approximation  problem  could  also  be  solved  by  using  the  Lp 
norm  with  p equal  to  a large  number.  Doing  so  may  result  in  sos* 
computational  advantages. 


5.  PHYSICAL  REALIZABILITY  OF  Z(s) 

The  impedance  function  Z(s)  derived  from  a TDR  voltage  v(t)  must 
be  physically  realizable  to  be  synthesized  with  resistors,  inductors, 
and  capacitors.  For  a Z(s)  given  by 


«*o  + 

to 

1 + 

a^s‘  + . . 

m 

. + a s 
m 

a 

, n 

bo  + 

b js  + 

b;'S‘  + . . 

. + b s 
n 

these  conditions  for  physical  realizability  can  be  expressed: 

a.  All  of  the  a^  and  b^  in  equation  (15)  are  real  and  positive. 

b.  The  poles  and  zeros  of  Z(s)  are  in  the  left  half  of  the 
s plane  or  on  the  imaginary  axis. 


c.  The  poles  and  zeros  on  the  imaginary  axis  are  simple. 

d.  The  real  part  of  is  not  negative  for  any  value  of  w. 


It  would  be  desirable  if  the  conditions  a through  d could  be 
guaranteed,  or  at  least  checked,  at  the  early  stages  of  the 
calculation.  However,  the  relationship  between  the  poles  and  residues 
of  the  SF.M  expansion  for  v(t)  and  the  poles  and  zeros  of  Z(s)  is 
complicated.  To  see  this  relationship,  consider  equation  (S)  for 
V(s)  : 


V(s) 


N A 

v ->  n ^ P(s) 

s - s Q(s)  ' 
n-1  n 


(16) 


3W.  H.  Chen,  Linear  Network  Design  and  Synthesis,  McGraw-Hill  Poo* 
Co.,  Inc.,  New  York  (1062). 
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where  P(s)  and  Q(s)  are  polynomials  of  degree  N-l  and  N,  respectively. 
Similarly,  the  SEM  expansion  of  v (t)  leads  to 


V> 


0„(s) 


(17) 


Substituting  equations  (16)  and  (17)  into  equation  (6)  yields 


2 (s) 


ZnQn (s) P (s)  _ P(s) 

P0<s)Q(s)  - Qo (s) P (s)  = Z°Q(s)  ' 


(IS) 


The  zeros  of  7.(s)  are  the  zeros  of  V(s)  plus  the  poles  of  Vu(s).  On 
the  other  hand,  the  poles  of  Z(s)  are  different  from  the  poles  and 
zeros  of  either  V(s)  or  Vj(s) . 


Usually,  the  generator  voltage  v()(t)  is  a unit  step  function  so 
that  V,t(s)  * 1/s.  Then  equation  (18)  can  be  simplified  to 


2 (s) 


2 O' 


sP(s) 
Q(s)  - 


sP(s) 


(19) 


This  equation  provides  an  easily  computed  relationship  between  the  SEM 
expansion  of  the  measurable  terminal  voltage  and  the  desired  impedance 
function.  Checks  for  the  satisfaction  of  conditions  a through  d have 
been  programmed  for  automatic  testing  with  the  digital  computer. 

A proce  lure  for  generating  a physically  realizable  impedance 
function  from  a given  V(s)  has  been  developed  and  tested.  The 
procedure  is  based  on  a theorem  in  Weinberg '*  and  is  described  in 
appendix  A.  The  results  of  numerical  tests  using  the  procedure, 
however,  have  not  been  encouraging.  In  general,  when  a nonrealizable 
impedance  function  is  obtained,  it  is  best  to  return  to  the  original 
v(t)  and  calculate  a new  SEM  approximation  with  different  values  of  N, 
K,  or  At. 

6.  SYNTHESIS  OF  LUMPED-ELEMENT  EQUIVALENT  CIRCUITS 

When  a physically  realizable  impedance  function  in  the  form  of 
equation  (15)  has  been  obtained,  a number  of  standard  circuit 
synthesis  procedures  may  be  employed  to  determine  an  RLC  network 
having  the  prescribed  input  impedance.  The  simplest  synthesis 


l*L.  Weinberg,  Network  Analysis  and  Synthesis , McCraw-Hill  Book  Co., 
Inc.,  New  York  (1962). 
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procedure  is  based  on  the  partial  fraction  expansion  of  the  impedance 
function  in  equation  (15).  This  procedure  has  failed  for  the  examples 
considered  because  the  individual  terms  of  the  partial  fraction 
expansion  are  not  all  physically  realizable.  Therefore,  more  general 
synthesis  procedures  such  as  those  of  Brune  or  Bott  and  Duffin  must  be 
used . 3 » 4 


The  Brune  procedure  has  been  selected  for  the  examples  that 
follow.  This  procedure  leads  to  the  least  number  of  circuit  elements 
for  a given  impedance  function.  A brief  description  of  the  Brune 
synthesis  procedure  is  given  in  appendix  B. 

To  illustrate  the  SFM  method  of  equivalent  circuit  synthesis,  a 
standard  gain  dipole  with  a \/4  balun  has  been  measured,  and  its 
equivalent  circuit  has  been  calculated.  The  original  data  obtained 
from  the  TDR  experiment  are  shown  in  figure  4.  An  SEM  approximation 
having  12  poles  is  shown  in  figure  5.  The  values  of  the  poles  and 
residues  are  qiven  in  table  II.*  The  maximum  deviation  of  this 
approximation  is  20  dB  below  the  0.5-V  signal  level,  and  the  total 
energy  in  the  error  siqnal  is  20  dB  below  the  energy  in  the  actual 
signal.  The  impedance  function  calculated  from  equation  (19)  is 


2 (s) 


ay  + a(S  + . . . + 

b o + b ) s + . . . + bUisl2 


(20) 


3t»\  H.  Chen,  Linear  Network  Design  and  Synthesis , McGraw-Hill  Book 
Co.,  Inc.,  New  York  (1962). 

**L.  Weinberg,  Network  Analysis  and  Synthesis , McGraw-Hill  Book  Co., 
Inc.,  New  York  (1962). 

* Si  nee  the  SEM  is  an  analytical  representation  of  the  experimental 
data,  the  poles  and  residues  qan  be  expressed  to  several  significant 
figures.  Maintaining  these  significant  figures  throughout  the 
calculations  insures  that  the  resulting  equivalent  circuit  will 
accurately  model  the  SEM  approximation  to  the  data.  Of  course,  the 
equivalent  circuit  model  for  the  antenna  cannot  be  more  accurate  than 
the  original  data. 


14 
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This  impedance  function  is  physically  realizable  as  can  be  seen  from 
the  coefficients  of  equation  (20),  the  locations  of  the  poles  and 
zeros  (table  III) , and  the  value  of  the  real  part  (fin.  6) . By  using 
Brune's  synthesis  procedure,  the  lumped-element  equivalent  circuit  in 
figure  7 is  obtained.  This  circuit  can  be  substituted  for  the  antenna 
in  any  network  analysis  code.  For  any  frequency  within  the  range  of 
validity  of  the  equivalent  circuit  (typically  several  decades),  the 
computed  voltaqes  and  currents  within  the  network  would  be  the  same  as 
if  the  antenna  were  attached.  For  receiver  applications,  a Thevenin 
representation  of  the  antenna  consisting  of  a source  and  the 
calculated  equivalent  impedance  circuit  can  be  used. 
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Figure  4.  Time-domain  ref lactometer  data  of  impedance 
of  dipole  with  1/4  balun. 


Figure  5.  Twelve-pole  singularity  expansion  method 

(HEM)  approximation  to  time-domain  ref lectometer 
(TDR)  data  of  dipole  antenna. 

TAPlf  II.  POLES  AND  RESIDUES  OF  SINGULARITY  EXPANSION  METHOD 
APPROXIMATION  TO  DIPOLE  DATA 
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TABLE  III 
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Figure  6.  Real  part  of  ^(s)  for  dipole. 

To  check  the  synthesis  procedure,  a circuit  that  represents  the 
TDR  experiment  has  been  analyzed  by  using  the  SPICE  computer  code. 
The  circuit  consisted  of  a 1-V  step  generator,  a 50-ohm  source 
resistor,  and  the  lumped-element  equivalent  circuit  (fig.  8) . The 
calculated  network  response  and  the  expected  response  (SEM 
approximation)  are  compared  in  figure  9.  The  comparison  is  excellent, 
and  the  equivalent  circuit  is  a very  good  representation  (errors  at 
least  20  dB  below  signal)  of  the  actual  antenna. 


.t 
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Figure  9.  Comparison  of  desired  (singularity  expansion 
method--SEM--approximation)  and  calculated 
responses  of  lumped-element  equivalent 
circuit  for  dipole. 

The  equivalent  circuit  for  a yaqi  antenna  has  also  been 
calculated.  Hie  TDR  data  and  the  10-pole  SEM  approximation  are  shown 
in  figure  10.  The  poles  and  residues  of  the  approximation  are  given 
in  table  IV.  Tbe  energy  of  the  error  signal  is  26  dB  below  the  energy 
of  the  actual  signal.  The  impedance  function  for  this  voltage  is 
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(21) 
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Fxamination  of  equation  (21),  the  data  in  table  V,  and  figure  11 
reveals  that  this  impedance  is  physically  realizable.  Tlte  Brune 
synthesis  procedure  yielded  the  equivalent  circuit  in  figure  12. 


Figure  10.  Ten-pole  singularity  expansion  method  (8F.M) 

approximation  to  time-domain  reflectometer  (TDK) 
data  of  yaqi  antenna. 

TABLE  IV.  POLES  AND  RESIDUES  OF  SINGULARITY  EXPANSION 
METHOD  APPROXIMATION  TO  YAG I DATA 
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l.umpt'il-o  l omen t equivalent  i ii'ui t tor  yaai  at\t  enna . 


A SPICE  analysis  of  the  TDK  response  of  this  lumped-element 
equivalent  circuit  yielded  the  result  shown  in  figure  13.  There  is 
excellent  agreement  between  the  expected  and  calculated  waveforms. 


Figure  13.  Comparison  of  desired  (singularity  expansion  method--SEM-- 
approximation)  and  calculated  responses  of  lumped-element 
equivalent  circuit  for  yagi. 

7.  TRANSFER  FUNCTION  CALCULATIONS 

The  terminal  properties  of  a linear,  time-invariant  antenna  can 
be  completely  described  by  a Thevenin  or  Norton  equivalent  circuit. 
These  circuits  consist  of  two  parts:  (1)  an  independent  source  that 

represents  the  open  circuit  voltage  or  short  circuit  current  of  the 
antenna  when  it  is  illuminated  by  an  electromagnetic  field  and  (2)  an 
equivalent  circuit  that  represents  the  antenna's  input  impedance  or 
admittance.  This  equivalent  impedance  or  admittance  circuit  can  be 
synthesized  by  the  techniques  described  above.  Furthermore,  the  SEM 
can  be  used  to  characterize  the  independent  source. 

The  SIX  characterization  of  the  antenna's  voltage  or  current  is 
obtained  by  calculating  the  impulse  response  or  the  frequency  domain 
transfer  function  of  the  antenna  in  the  form  of  equation  (7)  or  (8). 
The  characterization  is  calculated  by  finding  an  SEM  approximation  to 
the  measured  or  calculated  response  of  the  antenna  (open  circuit 
voltage  or  short  circuit  current)  to  a short  pulse  of  electromagnetic 
energy.2*5  This  SEM  approximation  is  an  analytical  estimate  of  the 

■M.  L.  Van  Blaricum  and  D.  H.  Schaubert,  An  Experimental  Transient 
Transfer  Function  via  Prony's  Method,  Proc.  1976  USNC/URSI  Meeting, 
Amherst,  MA  (October  1976) , 97. 

5C.  H . Schaubert , Measurement  of  the  Impulse  Response  of 
Communication  Antennas , Harry  Diamond  Laboratories  TR-1832  (November 
1977) . 


antenna's  impulse  response,  and  it  may  be  convolved  with  any  incident 
waveform  to  obtain  the  resulting  terminal  voltage  or  current. 
Therefore,  the  SEM  can  be  used  to  generate  equivalent  circuits  that 
completely  characterize  the  terminal  properties  of  an  antenna. 

8.  CONCLUSIONS 

The  singularity  expansion  method  (SEM)  has  been  shown  to  be  an 
effective  tool  for  synthesizing  lumped-element  equivalent  circuits 
from  experimentally  derived  time-domain  reflectometer  (TDK)  data.  The 
poles  of  the  SEM  expansion  are  calculated  by  using  the  standard  Prony 
algorithm.  The  residues,  however,  are  calculated  by  an  iterative 
scheme  that  yields  a uniform  norm  approximation  instead  of  a least 
squared  error  norm  approximation.  The  uniform  norm  provides  much 
better  control  of  the  early  time  errors  and  leads  to  a better  overall 
approximation.  The  increased  computation  time  required  for  the 
uniform  norm  approximation  is  readily  justified  because  this  time  is 
not  large  compared  to  the  total  time  for  calculating  the  equivalent 
circuit.  Furthermore,  the  time  spent  in  calculating  the  SEM 
approximation  is  repaid  by  the  final  circuit,  which  produces  exactly 
the  voltaqe  response  prescribed  by  the  SEW  approximation. 

Once  the  physical  realizability  of  the  calculated  impedance 
function  has  been  verified,  a standard  circuit  synthesis  procedure, 
such  as  that  of  Rrune  or  Pott  and  Duff in,  may  be  used  to  obtain  the 
desired  RLC  equivalent  circuit.  In  general,  the  partial  fraction 
expansion  of  F(s)  will  not  yield  subcircuits  that  are  physically 
realizable. 

When  combined  with  the  previously  developed  technique  for 
measuring  the  impulse  response  of  antennas,  SEM  synthesis  of 
equivalent  impedance  circuits  is  a powerful  and  generally  applicable 
tool  for  converting  experimentally  obtained  data  into  Thevenin  and 
Norton  equivalent  circuits  for  use  in  wide-bandwidth  analytical 
modelling  of  transmitter  and  receiver  systems. 
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APPENDIX  A.— PHYSICAL  REALIZABILITY  VIA  MULTIPLICATIVE  TACTOR 


Theorem  6.S  in  Weinberg1  states  that  7(s)  ■ 7.  , 1 1 + R(s))/|l 
R t s')  1 is  a physically  realizable  function  if  and  only  if 

1.  R(s)  is  a real  rational  function  with  no  poles  in  the  riqht 
half  plane  or  on  the  imaginary  axis,  and 

2.  |R(1  .0  | _v  1 for  all  real  u>. 

This  theorem  and  equation  (18)  in  the  main  body  of  this  report 
can  I^e  used  to  obtain  an  impedance  function  7^(s)  that  depends  on  the 
parameter  i,  0 < a s 1.  For  a *•  1 , 7 (s)  equals  the  impedance 
correspondinq  to  R(s)  . For  0 «:  a < 1,"  Za(s)  is  the  impedance 
correspondinq  to  a reflection  coefficient  aR(s) . 

The  proposed  scheme  for  calculatinq  a physically  realizable 
im|iedance  function  requires  that  the  singularity  expansion  method 
(SEM)  approximation  to  V(s)  and  V0(s)  be  manipulated  to  obtain 
R(s)  « N(s)/D(s),  where  N(s)  and  D(s)  are  polynomials.  ihen 


r* 

I* 


(s) 


1 + aR(s) 
l’  1 - UR(s) 


(A-l) 


can  be  calculated  for  various  a,  and  the  larqest  a < 1 for  which  7 (s) 
is  physically  realizable  is  the  "best"  impedance  function.  ("Best"  is 
used  in  the  sense  that  the  poles  and  zeros  of  the  reflection 
coefficient  are  those  obtained  from  V(s)  and  v,i(s>)  and 
;k  (j,)  : - > R(j.)  < 1.) 

The  impedance  function  Za(s)  can  lie  related  directly  to  the  SEM 
expressions  for  V(s)  and  Vl,(s).  Equation  (?)  in  the  main  body  of  the 
report  implies  that 
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l.  Weinberg , Network  Analysis  and  Synthesis , McGraw-Hill  nook  Co., 
Inc.,  New  York  (1**2) . 


(1  - a)Q(s)P(,  (a)  + 2ciP(s)Q(l  (s) 


VS)  Z°  ( 1 + u ) Q ( s ) P o ( s ) - 2rtP(s)Q0(s)  ■ (A_4: 

For  step  function  excitation,  V0(s)  ■ i/s,  and  equation  (A-4)  reduces 
to 

(1  - u)Q(s)  + 2otsP(s) 

Z (s)  = Z0— 

* U(1  + u)Q(s)  - 2asP(s)  * 


APPENDIX  B. --BRUNE  SYNTHESIS  OF  RLC  NETVORKS 


A detailed  description  of  the  Brune  synthesis  procedure  can  be 
found  in  many  texts.  1,2  The  followinq  is  a brief  summary  of  the 
procedure . 

Tlie  Brune  synthesis  procedure  is  based  on  the  tee  representation 
of  a transformer  (fig.  B-l) . In  the  tee  representation,  one  of  the 
inductances  L or  Lc  is  negative,  but  the  network  is  still  considered 
physically  realizable. 


1»  L« 

__rrm_^_rrrn 


Ip  - L»+  Li 
M - lk 
Lt  " L»+  Lt 

Figure  B-l.  Equivalent  representations  of  transformer. 

Starting  with  the  nth  order  impedance  function, 

n n- 1 , , 

as  + a s + . . . + a() 

Z(a>  = -S- -n"J-TiV- . 

b sn  + b s + . . . + bo 

n n- 1 

a single  cycle  of  the  synthesis  procedure  yields  the  network  in 
figure  B-2,  where  Z i,(s)  is  of  order  n-2.  The  Brune  cycle  for  reducing 
Z(s)  consists  of  four  steps: 


l£.  Weinberg,  Network  Analysis  and  Synthesis , McGraw-Hill  Book  Co., 
Inc.,  New  York  (1962). 

' w . H.  Chen,  Linear  Network  Design  and  Synthesis,  McGraw-Hill  Book 
Co.,  Inc.,  New  York  (1962). 
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APPENDIX  B 

1.  Remove  Rj  «=  m*n  Re  Z(jou)  (fig.  B-3)  . 

2.  Remove  Lj  = im  ZfjuijJ/Uj. 

3.  Remove  the  shunt  arm  L C , which  produces  zeros  in  Z0 (s) 
at  s = ijiii). 

4.  Remove  L3,  which  causes  Z3<s)  to  have  a pole  at  infinity. 


Figure  B-2.  Realization  of  single  Brune  cycle. 
Hi  Z(M 


Figure  B-3.  Typical  resistance  function  for  Brune 
synthesis  procedure. 
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appendix  b 


•rtiese  four  steps  require  polynomial  manipulations  that  can  be  readily 
proq rammed  on  the  digital  computer.  The  equivalence  in  figure  B-l  may 
be  used  to  convert  a Brune  circuit  containing  a negative  inductance 
into  a transformer  as  shown  in  figure  B-4. 


Figure  B-4. 


Final  Brune  circuit  with  transformer. 
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